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ABSTRACT 

We study the role of the guide field in relativistic magnetic reconnection in 
a Harris current sheet of pair (e ) plasmas, using linear theories and particle- 
in-cell (PIC) simulations. Two-dimensional PIC simulations exhibit the guide 
field dependence to the linear instabilities; the tearing or reconnection modes 
are relatively insensitive, while the relativistic drift-kink instability (RDKI), the 
fastest mode in a relativistic current sheet, is stabilized by the guide field. Par- 
ticle acceleration in the nonlinear stage is also investigated. A three-dimensional 
PIC simulation demonstrates that the current sheet is unstable to the RDKI, 
although a small reconnection occurs in the deformed current sheet. Another 
three-dimensional PIC simulation with a guide field demonstrates a completely 
different scenario. Secondary magnetic reconnection is triggered by nonlinear 
coupling of oblique instabilities, which we call the relativistic drift-sausage tear- 
ing instability. Therefore, particle acceleration by relativistic guide field recon- 
nection occurs in three-dimensional configuration. Based on the plasma theories, 
we discuss an important role of the guide field: to enable non-thermal particle 
acceleration by magnetic reconnection. 

Subject headings: acceleration of particles — magnetic fields — plasmas — in- 
stabilities — relativity 

1. INTRODUCTION 



Magnetic reconnection in collisionless plasmas is of strong interest in high-energy astro- 
physical places such as active galactic nuclei (AGNs; di Matteo (1998); Birk et al. (2001)), 
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extragalactic jets (Romanova & Lovelace 1992; Lesch & Birk 1998 Larrabee et al. 2003), 



pulsar winds (Michel||1982 


1994; Coroniti||1990 Lyubarsky & Kirk||2001 Kirk & Skjseraasen 


2003 


), gamma-ray bursts 


Drenkhahn 2002 


Drenkhahn & Spruit 2002 Uzdensky & Mac- 



Fadyen 


2006 


) and, soft gamma repeaters ( 


Thompson & Duncan 


1995 


2001 


Lyutikov 


2003 



2006). Since it rapidly releases stored magnetic energy into plasma kinetic energy, recon- 
nection is considered as a possible underlying mechanism, to explain particle acceleration or 
bursty emission signatures in these sites, regardless of plasma composition (pair plasmas or 
ion-electron plasmas). 

For instance, in the case of the Crab pulsar, magnetic reconnection is considered in the 
relativistic radial outflow of pair plasmas. Since the central neutron star is a fast oblique 
rotator, its strong magnetic fields (~ 10 12 G) are highly striped so that field lines are almost 
toroidal inside the flow. Magnetic reconnection is expected to occur in field reversal config- 



uration of toroidal fields, in order to dissipate the magnetic energy (Coroniti 1990). In the 
AGN context, magnetic reconnection is expected in association with magnetic loops from 
a differentially rotating acceleration disk. In this situation, the field configuration involves 
magnetic fields perpendicular to the antiparallel component. Indeed, reconnection in sheared 



configuration is proposed as an acceleration site of MeV/GeV electrons (Lesch & Birk 1997 



Schopper et al. 1998 Nodes et al. 2003), due to the field-aligned electric field En. In the 
case of soft gamma repeaters, a giant flare is expected in ultrastrong magnetic environment 
(~ 10 15 G) around central neutron stars (magnetars). Although there are various theoretical 
models, it is expected that giant flares involve magnetic reconnection in the relativistic pair 
plasma environment. 

However, the fundamental mechanism of relativistic magnetic reconnection, as well as 
the conventional non- relativistic counterpart, is far from being understood. From the mag- 
netohydrodynamics (MHD) viewpoint, steady state reconnection models have been extended 



to the relativistic regime (Blackman & Field 1994; Lyutikov & Uzdensky 2003 Lyubarsky 



2005), although they do not fully agree. The relativistic resistive MHD simulation demon- 
strated mildly-relativistic Petsheck reconnection ( |Watanabe fc Yokoyama 2006). From the 
kinetic viewpoint, series of particle- in-cell (PIC) simulations demonstrated that relativistic 



magnetic reconnection is a powerful acceleration engine (Zenitani & Hoshino 2001, 2005b 



2007 Jaroschek et al. 2004b), primarily due to the direct particle acceleration by the recon- 
nection electric field, which is perpendicular to the reconnecting magnetic field lines. In fact, 
obtained plasma energy distribution is highly non-thermal, approximated by the power-law 



distribution with an index of —1 ~ —3 ( 


Romanova & Lovelace 


1992 


Zenitani & Hoshino 


2001 


Larrabee et al.||2003 


Jaroschek et al.||2004b 


). Thus, magnetic reconnection seems to be 



a favorable source of non-thermal particles at the sites of synchrotron radiation (Jaroschek 



et al. 2004a). However, it was reported that the current sheet configuration is unstable to 
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the relativistic drift kink instability (hereafter the RDKI), which occurs in the plane per- 



pendicular to the reconnection plane (Zenitani & Hoshino 2005a). The RDKI grows more 



rapidly than reconnection, and it mainly contributes to plasma heating ( Zenitani & Hoshino 



2007). Thus, plasma heating by the RDKI overwrites the reconnection scenario; nonthermal 
particle acceleration is unlikely to occur in antiparallel field configuration in relativistically 
hot pair plasmas. 

Pair plasma reconnection has also attracted recent attention as an equal-mass-limit ex- 
ample of ion-electron reconnection. For example, there is a long-standing problem regarding 
the reason collisionless reconnection occurs "faster" than predicted by the MHD theories. It 
has been widely argued that Hall physics plays an essential role in maintaining fast recon- 



nection (the GEM reconnection challenge; Birn et al. (2001); Shay et al. (2001)), but fast 



reconnection in non-relativistic pair plasma shows a counter-evidence because it does not 



involve Hall effects (Bessho & Bhattacharjee 2005). Instead, although its physical interpre- 



tation needs refinement, it seems that the off-diagonal part of the pressure tensor (Hesse et 



al. 1999) accounts for the electric field for fast reconnection (Bessho & Bhattacharjee 2007 



Hesse & Zenitani 2007). Furthermore, it has recently been argued that late-time dynamical 



structure may regulate fast reconnection (Daughton & Karimabadi 2007) 



Guide field configuration, which contains a perpendicular magnetic field to the antipar- 
allel components, is an important generalization for studying the reconnection problem in 
a shear or twisted configuration. Reconnection with the uniform guide field has long been 



studied for better understanding of the reconnection structure (Drake & Lee|1977 


Katanuma 


& Kamimura 1980 


Hoshino 1987 


Horiuchi & Sato 1997 Hesse et al. 1999 2004) and for 


applications to solar flares and to the Earth's magnetopause (Shibata et al. 


1995 Quest 



& Coroniti 1981 Sonnerup et al. 1981). Recent three-dimensional (3D) PIC simulations 



(Scholer et al. 2003 Drake et al. 2003 Ricci et al. 2003; Silin k Buchner 2003 Pritchett 



& Coroniti 2004) discussed reduction of the cross-field activities (Scholer et al. 2003 Silin 



& Buchner 2003), detailed reconnection region structure (Pritchett & Coroniti 2004), and 



enhanced electron acceleration (Drake et al. 2003 Ricci et al. 2003). We note that the 



acceleration mechanism differs from the one in the antiparallel counterpart because the re- 
connection electric field is no longer perpendicular to magnetic lines; guide field reconnection 



involves a parallel electric field E\\ to accelerate electrons (Ricci et al. 2003). In relativistic 



pair plasmas, Zenitani & Hoshino (2005b) presented the idea that the guide field stabilizes 
the RDKI and that the new oblique instability triggers secondary magnetic reconnection. 
Nonthermal particle acceleration by magnetic reconnection seems to be a likely process in a 
relativistic current sheet once again, under the guide field condition. 



The purpose of this paper is to comprehensively discuss the role of the guide field 



-4- 



in relativistic pair plasma reconnection. Extending our recent papers (Zenitani & Hoshino 



(2005b) and Zenitani & Hoshino (2007), hereafter referred as ZH05b and ZH07, respectively), 



we investigate the linear and non-linear development of a relativistic current sheet in pair 
plasmas, by using fully self-consistent kinetic PIC simulations. The conclusion of this paper 
is that relativistic guide field reconnection is a favorable process for nonthermal particle 
acceleration, from the viewpoint of detailed plasma theory. 

This paper consists of the following sections. In §2 we describe the simulation setup. 
In §3 we study the two-dimensional evolution of the current sheet with the guide field. We 
investigate how the guide field affects the linear instabilities in §3.1, and then we study how 
particles are accelerated in nonlinear guide field reconnection in §3.2. In §4, we study the 
3D evolution of an anti-parallel current sheet. In §5, we study the 3D evolution of a current 
sheet with a guide field. We introduce new linear instabilities in the oblique directions in 
§5.1, and then we discuss the nonlinear development of the current sheet in §5.2, in relevance 
to the two-dimensional counterparts. Section 6 contains discussion and the summary. 



2. SIMULATION MODEL 

Throughout this study, we employ a relativistic Harris model as an initial configuration. 
Magnetic field and plasma distribution functions are set in the following way: 

B = B tanh(z/X)x + aB y, (1) 

n cosily/A) - ls {e - I3 s mcu y ) _ j_ 

Js A-nm*cTK 2 {mc 2 /T) P[ T J + 4irm 2 cT bg K 2 (mc 2 /T bg ) P [ T bg h K ) 

where B is the antiparallel magnetic fields, a is the relative amplitude of the guide field, n 
is the plasma number density of the current sheet in the proper frame, T is the plasma tem- 
perature including the Boltzmann constant in the proper frame, m is the positron/electron 
mass, c is the light speed, K 2 (x) is the modified Bessel function of the second kind, the sub- 
script s denotes the species (p for positrons and e for electrons), c[3 s is the drift speed of the 
species, and u is the four velocity. The parameters T and hq are defined in the proper frame. 
T bg and n bg are the temperature and the number density of background plasmas. In this 
paper, we choose (3 P = 0.3, (3 e = —0.3 (relevant Lorentz factor is 7/3 = 1.048), T/mc 2 = 1, 
T bg /T = 0.1 and 77,^/(79710) = 0.05. Notice that the Harris model with the uniform guide 
field exactly satisfies an equilibrium. We assume that B y is negative, but we will observe 
the same (mirror) results in the positive cases, because of the mathematical symmetry in 
pair plasmas. We use 3D particle-in-cell (PIC) code with periodic boundaries in all three 
directions. Because of the field reversal, we set double current sheets in the z direction. The 
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typical scale of the current sheet A is set to 10 grids. The size of the simulation box depends 
on the problem. Their sizes are presented in Table [TJ These conditions are same as our 
recent studies; runs R3 and D3 are identical to ones in ZH07 and runs 3D-A and 3D-B are 
identical to run A and run B in ZH05b. 

In many cases, we impose small artificial electric fields around (x,y, z) = (0, 0, ±3A) in 
the very early stage of simulation runs, in order to "force" reconnection around the center 
of the main simulation domain. Typical spatial ranges for the artificial fields are set to 
(Ax, Az) ~ (±2A, ±A) in two-dimensional cases and (Ax, Ay, Az) ~ (±2A, ±2A, ±A) in 3D 
cases. This force is often strong enough and then reconnection process immediately breaks 
up soon after short linear growth stage. In Table [TJ the letter F means that we observe 
such sudden breakup of forced reconnection. S* means that we observe quasi-spontaneous 
evolution, because the trigger force is not strong enough. In these cases, it takes some time 
until nonlinear reconnection breaks up, and we observe the linear growth of the instabilities 
for a while. In S cases, we set no trigger force. In two-dimensional runs for the RDKI, we 
do not need the trigger force because the instability quickly grows from thermal noise. 



3. GUIDE FIELD EFFECT TO THE TWO-DIMENSIONAL PROCESSES 

3.1. Linear Growth Rates 

First, we investigate the guide field effect on the two-dimensional reconnection insta- 
bilities: relativistic reconnection in the x-z plane and the RDKI in the y-z plane. We 
carry out a series of two-dimensional simulations: five cases of reconnection (a = B v /Bq = 
0, -0.5, -1, -1.5, -5) and four cases of the RDKI (B y /B = 0, -0.25, -0.5, -1). Then, we 
calculate the linear growth rates of the fastest growing modes from the perturbed magnetic 
fields 5B in the neutral plane. In reconnection cases, it is sometimes difficult to pick up 
the "linear" mode, because the linear stage is rather short. It seems that the linear stage 
takes longer time, when the initial artificial impact is zero or weak. Therefore, for the most 
difficult case of B y /B = —1.5 (run R3-C), we also carry out the relevant spontaneous run 
(run R3-D) to confirm the linear growth rate. We notice that the linear growth rate is not 
identical to the energy release rates in the nonlinear stage. Once the nonlinear reconnection 
breaks up, the linear perturbations quickly cascade into the longer-wavelength perturbations, 
which grows more explosively, 2-3 times faster than the linear or quasi-linear modes. 

Along with the PIC simulations, we evaluate the linear growth rates of the two-dimensional 
instabilities by using eigenvalue analysis, in the same way as our previous work (Appendix 
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V S , 9Ps 

-^{7 s q s n s E ■ v s H 



(3) 



B in ZH07). In this analysis, we use the relativistic fluid equations 
7 2 d ( 

-^{p s + e s )(^T + v, ■ V)u fl = - Vp s + "f s q s n s \E + — x ti J - — {j s q s n s ±J ■ v s + 

where p is the isotropic plasma pressure and e is the internal energy that contains the 
rest mass energy. In addition, we use the continuity, Maxwell equations, and adiabatic 
gas condition to close equations. We assume the perturbation for arbitrary wavevector 
k = (k x ,k y ); Sf oc 5f(z)exp(ik x x + ik y y + cjjt), where Ui is the linear growth rate (the 
imaginary part of the wave frequency), and then we solve the linearized equations as an 
eigenvalue problem. 

Let us discuss two-fluid theory in the x-z plane (k y = or d/dy = 0). By ignoring the 
relativistic drift speed effect (third term of the right-hand side of eq. [3]), simple assumptions 
(Qp = ~Qe = Q, 1p = 7e = 7, P P = Pe = P, e p = e e = e, and n p = n e = n) yield the y 
component of the generalized Ohm's law 

7(p + e 



V 

E+ — x B 

c 



2 c 2 qn 



d_ 

dt Vpy 



d_ 

dt Vey 



(4) 



where V = (1/2) (v p + v e ) is the bulk velocity of pair-plasma fluids. This indicates that the 
fluid inertia of electrons and positrons account for the effective resistivity to break down the 
frozen-in condition and to drive the tearing instability. The Hall terms are not included, 
because they vanish in electron-positron plasmas. As discussed in the previous works on the 



relativistic tearing instability (Zelenyi & Krasnoselskikh 1979) and the conventional tearing 
instability, we obtain purely-growing modes in the wavevector range of < k x X < 1. We 
employ the maximum growth rates in the above parameter range as theoretical eigen growth 
rates of the reconnection mode. In the y-z plane (k x = or d/dx = 0), Daughton (1999) 
demonstrated that two-fluid theory and Vlasov theory are in good agreement to explain 
the drift kink instability in ion-electron plasmas, in long- wavelength range of k y X < 0.7. 
Similarly, we know that our two-fluid theory is in agreement with PIC simulations in long- 



wavelength rage of the RDKI; k y X < 0.7 (Zenitani & Hoshino 2005a, 2007). Therefore, we 
employ k y X = 0.7 to estimate the typical eigen growth rates of the RDKI. We keep in mind 
that our theory assumes isotropic fluids and that it does not contain any kinetic effects nor 
wave-particle interactions. These approximations will break down in a thin current sheet, 
in which gyro radii is much larger, and so kinetic-based theories are more favorable for 



describing instabilities (e.g. Vlasov analysis by Brittnacher et al. (1995)). However, in the 
scope of the present paper, two-fluid theory gives us a plausible estimate of the growth rates. 

The obtained growth rates of the reconnection mode are presented in black lines in 
Figure [I] as a function of \B y /B \. In PIC simulations, the fastest reconnection mode (Fig. 
[Tj black line) has wavenumbers of k x X ~ 0.4—0.7, which correspond to those of the relativistic 
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tearing instability. The maximum eigen growth rates of the relativistic tearing instability 
are presented with the dashed line in Figure [T] Both growth rates are in agreement. It 
seems that the growth rates are relatively insensitive to the guide field. In the extreme guide 
field case (\B v /Bq\ = 5.0), the growth rate is still the half of that of the antiparallel case. 
This tendency is consistent with reconnection studies in ion-electron plasmas; an analysis 
(Brittnacher et al. 1995) and PIC simulations (Horiuchi & Sato 1997) reported that the 
linear growth rate decreases by a factor of 2 or 3 under strong guide field conditions (up to 
\B y /B \ =4). It was also reported that the guide field reduces the quasi-steady reconnection 



rate by a factor of ~ 2 in PIC simulations and Hall MHD simulations (Pritchett & Coroniti 



2004 Huba 2005). 



The gray lines in Figure [T] show the growth rates of the RDKI (gray line by PIC simu- 
lation; gray dashed line by the theory). The RDKI is rather sensitive to the guide field. In 
PIC simulations, the RDKI spontaneously arises from the thermal noise before t/r c = 100 
under no/weak guide field conditions (Runs D3 and D3-A; \B y /B \ = 0,0.25). However, we 
could not observe any sign of the RDKI before t/r c = 200, when the guide field is strong 
enough, \B y /B \ > 0.5. We also failed to obtain the unstable eigen solutions when guide 
field amplitude exceeds a threshold; \B y /B \ > (0.45 — 0.5). In the physical sense, since the 
wavevector is parallel to the guide field lines, the magnetic tension of the guide field prevents 
the RDKI from growing. Previous ion-electron work (Pritchett & Coroniti 2004) reported 
that the guide field B y /B = 0.5 is sufficient to stabilize the kink modes. 

In summary, in two-dimensional linear regime, the RDKI grows faster than the recon- 
nection or the relativistic tearing instability under no/weak guide field conditions. However, 
the situation changes when we introduce a finite amount of the guide field. The reconnection 
is rather insensitive to the guide field, but the RDKI is suddenly stabilized by the guide field. 
Therefore, reconnection will dominate in the guide field conditions of \B y /B Q \ > 0.5. 



3.2. 2D Guide Field Reconnection 

Next, we study the nonlinear development of the relativistic guide field reconnection. 
Among several simulation runs, we present the results of run R3-C (\B y /B \ = 1.5), in which 
the guide field features are well developed. 

We force magnetic reconnection around the center and then reconnection quickly breaks 
up. Figure [2] shows snapshots at t/r c = 100, where r c = A/c is the light transit time. Figure 
[2]a shows plasma density, average flow and magnetic field lines. The maximum outflow speed 
is ~ 0.37c along the x-axis and it is substantially slower than the antiparallel case of ~ (0.7- 
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0.8)c. The composition of plasma flow depends on the quadrants in the reconnection region; 
electron flow into the (+x, —y) direction in the (+x, +z) quadrants, positron flow into the 
(—x, +y) direction in the (— x, +z) quadrants, electron flow into the (— x, —y) direction in the 
(+x, —z) quadrants and positron flow into the (+x, +y) direction in the (+x, —z) quadrants 
are dominant. The X-type layer along the separatrix region is also characteristic. Figure 
[2^ shows the ^-component of the positron current. The current is strong around the O-type 
region, while it is weak or reversed in the very center of the O-type region, where the plasma 
density is at a maximum. Around the reconnecting region, there is a broadened structure of 
weak y-current. The most characteristic thing is the inclined current layer along the lower 
side of the separatrix, which is indicated by the dashed lines. This inclined layer corresponds 
to one of the X-type density regions in Figure [2]a. The electron current structure is upside 
down; it is slightly inclined in the counter-clockwise direction and it is on the upper side 
of the separatrix. In ion-electron plasmas, several simulations of guide field reconnection 



reported an inclined electron current layer (Horiuchi & Sato 1997 Drake et al. 2003; Hesse 



et al. 2004 Pritchett & Coroniti 2004) near the X-type region. However, the ion current 
layer is hardly recognized due to the ion's large spatial scale. In our simulation, we observe 
both current layers for positrons and electrons. Positrons are driven by the reconnection 
electric field E y around the X-type region, and then they tend to escape in the lower right 
(or upper left) direction along 3D reconnected field lines. Similarly, electrons tend to escape 
from the X-type region; they tend to escape in the upper right (or lower left) direction. 
Therefore, there is charge non- neutrality around the outflow region (Fig. |2]c). Positrons are 
localized in the lower right side of the reconnection region. Figure [2]c? shows the reconnection 
electric field E y . Its typical amplitude is E y /B ~ (0.05 — 0.1) around the X-type region and 
it is substantially weaker than antiparallel case; E v /Bq ~ (0.2 — 0.3). The E x component is 
comparable or slightly larger than E y , but consistent with plasma inflows. Its typical value 
is E X /B Q ~ 0.1 in the upper inflow region and E X /B Q ~ —0.1 in the lower inflow region. The 
vertical electric field is substantially large: E z /Bq ~ 0.7 in Figure [2]e. It is consistent with 
the outflow jet and guide magnetic field B y , which is compressed around the O-type region. 
The non-neutral charge distribution is consistent with the E z structure, too. 

Figure [3] shows the energy spectra over the main simulation domain. As reconnection 
evolves, particle acceleration takes place and then the high-energy tail continues to grow. 
The late-time spectrum at t/r c = 200 is quasi-stable. Although it is difficult to discuss the 
spectral index in such a temporally and spatially limited system, the power-law index is 
~ —2.9 in the last stage. 

In order to study particle acceleration, we track the trajectories of the highest energy 
positrons, whose energy exceeds 40mc 2 at t/r c = 100 (indicated by gray shading in Fig. 
[3]). Over the entire simulation domain they are only found along the inclined current layer, 



-9- 



both in the upper left side and in the lower right side in the reconnection region. The 
high-energy electrons are found in the other current layer. Gray points in Figure |4]a show 
spatial distribution of selected positrons around the lower right current layer. We pick up 
two typical trajectories from t/r c = to t/r c = 200, which are indicated by solid lines in 
Figure |4} We call them positron A and positron B, respectively. The diamond or circle signs 
show their positions at t/r c = 100. Positron A belongs to the majority. Starting from the left 
outside the presented region, the positron visits the A-type region and then escapes toward 
the lower right, gyrating around the reconnected field lines. The path is along the inclined 
positron current layer. While it travels a long distance in the ^/-direction, the particle mainly 
gains its energy from the reconnection electric field E y near the A-type region. Notice that 
the scaling of the ^-direction is extremely larger than the x- and z-directions in Figure |4]&. 
The energy gained from E x and E z is substantially low. Positron B is a rare sample. It 
starts near the center and then it moves into the A-type region. As reconnection goes on, it 
travels in the ^-direction along the A line in the x-z plane because it is trapped by the guide 
field By around the A line. At tjr c = 200, it still rotates around the A line and its energy 
increases from 40.8mc 2 (t/r c = 100) to 75.5mc 2 (£/r c = 200). The acceleration continues until 
it finds a way to escape toward the upper left or the lower right. 

In the antiparallel reconnection, the reconnection electric field E y is a main player of 
particle acceleration, too. Note that the reconnection electric field is perpendicular to the 
magnetic field lines in the antiparallel case, EB ~ 0. Therefore, particle acceleration occurs 



in the acceleration region where the frozen-in condition is never satisfied (\E\ > \B\; Zenitani 



fc Hoshino|fl200ll )), or in the piled-up regions where the electro-magnetic fields have strong 



peaks (ZH07). In the present case of guide field reconnection, particle acceleration occurs 
because E ■ B ^ around the reconnection region. The guide field B y traps the particles, 
preventing them from escaping in the x- or z-directions, so that they are accelerated by 
E y . The particle acceleration occurs even when they gyrate around the magnetic field lines, 
because it relies on the parallel electric field E\\. 

We have also compared the nonlinear evolutions of two-dimensional reconnection runs 
with various guide field amplitudes. Generally speaking, as the guide field increases, sig- 
natures of guide field reconnection become more apparent; electron/positron flows become 
separated, and then the inclined current layers and charge separation structure appear. The 
reconnection electric field also becomes weaker because outflow jets becomes slower. In ad- 
dition, in spontaneous or nearly spontaneous cases, it takes a longer time until the nonlinear 
reconnection process breaks up, as the stronger guide field is imposed. 



Figure 15 compares the acceleration efficiency in various reconnection runs. The ratio of 
the nonthermal kinetic energy to the total kinetic energy are presented as a function of time. 
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The ratio is calculated from the energy spectra in the same way as ZH07 (Appendix A in 
ZH07); we assume an equivalent energy spectra of the relativistic Maxwellian distribution, 
which carries the same amount of kinetic energy as the energy spectra in simulations, and 
then we obtain the nonthermal kinetic energy by subtracting the equivalent thermal spectra 
from the simulation spectra in their high-energy tail. This method may underestimate 
the amount of the nonthermal energy, but we can quantitatively compare the acceleration 



efficiency in a simple way. In Figure 15, the simulation time for run R3-A (\B y /B \ = 0.5), in 
which the trigger force is slightly weak, is re-arranged by At = 35r c so that we can directly 
compare it to the other cases. Here we focus on the two-dimensional (2D) reconnection 
runs; we discuss the 3D results later again in §5.2. Roughly speaking, the acceleration 
efficiency seems to decrease as the guide field increases. The antiparallel run (run R3) is the 
most efficient accelerator among the three runs. The maximum accelerated energy is also 
a function of the guide field during the similar timescale (See Table [T|. This is due to its 
strong reconnection electric field E y ~ 0.3Bq, and the reconnection electric field E y becomes 
weaker in the guide field cases. In run R3-C, the "acceleration region" nearly smears out 
due to the weak reconnection field E y ~ 0.1-B . On the other hand, the trapping effect of 
the guide field contributes to the acceleration efficiency. Even though particle acceleration 
by E± becomes less active, parallel acceleration by E\\ works instead. Thus, guide field 
reconnection (run R3-C) still produces a substantial amount of nonthermal energy ~ 15%. 
Considering the limitation of our simple method, more kinetic energy will be carried by 
nonthermal particles. In general, guide field reconnection seems to be an efficient particle 
accelerator, which releases more than 15%-30% of the plasma energy into the nonthermal 
kinetic energy. 

In run R3-A, we notice that the secondary particle acceleration takes place in the very 
late stage (at t/r c > 250), and then particles are accelerated up to 7 ~ 200. The spectral 
index evolves to a harder value of ~ —2.2 in the very late stage (at t/r c = 400). Since 
periodic boundary effects arise in that late stage, particle acceleration may be enhanced due 
to interactions with the multiple reconnection structures. The secondary acceleration is only 
found in the weak guide field case (run R3-A), and so its long-term evolution needs further 
investigation. 

When we compare runs R3-D and R3-C (\B y /Bo\ = 1.5), the system evolution in run 
R3-D delays by approximately 145t c -150t c in terms of the global energy distribution. In 
run R3-D, we observe six tearing islands around t/r c = 200, and then they start to collide 
with each other until one reconnection region dominates around t/r c = 300. The O-type 
islands often have onion-ring density structure, and the A-type current layers connect to 
the outermost rings. However, the global energy evolution and the acceleration efficiency 
are almost the same. They seem to be rather insensitive to the initial evolution whether 
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the system evolves from forced single reconnection or from multiple tearing islands, in the 
moderate guide field cases of \B y /B \ ~ 1.5. 

4. 3D EVOLUTION WITH NO GUIDE FIELD 

In this section, we study the 3D evolution of the current sheet in antiparallel configu- 
ration. We look at the simulation results of run 3D-A. As stated in §2, we imposed small 
external electric fields around the center, but the trigger fields are not strong enough to force 
reconnection. 

We briefly review the linear evolution, which was discussed in our previous paper 
(ZH05b). Figure |5]a is a snapshot at t/r c = 80. The current sheet is between the two 
gray surfaces and the plasma density at the neutral plane (z = 0) is projected into the 
bottom wall, with color shading from black (empty) to red (dense; n = 1.2n ). Figure [6] 
presents the E y structure at t/r c = 80. The red regions have the positive polarity of E y > 0, 
while the blue regions are negative: E y < 0. Apparently these profiles exhibit a quasi 2D 
evolution of the RDKI. One can see the polarity change along the x-axis in the x-z plane 
(the right surfaces in Fig. [6]), but this is due to small y-displacement of the structure. The 
wave-number of the RDKI is k y X ~ 0.74 (mode 3), while mode 4 is observed in the relevant 
2D run D3. However, both mode 3 and mode 4 are reasonable with the linear theories. The 
observed growth rate T c Ui = 0.06 is slightly slower than the expected rate T c Ui ~ 0.1, but it 
is still faster than that of the relativistic tearing instability. 

Figures [5]6 and|5]c show the nonlinear development of the current sheets at t/r c = 110 
and 140. The current sheet is strongly folded at t/r c = 110, and then its wave fronts start to 
collide each other. At t/r c = 110, the gray density surfaces are reset to n = l/3no, so that 
we can see the low-density hole around the center. We discover that magnetic reconnection 
takes place in this central hole. Figure [5](i is a zoomed-in view around the neutral plane at 
t/r c = 110. The gray lines are magnetic field lines, traced from the following six start points: 
(x,y,z) = (— 5, 0, — 0.5), (— 5, 2, — 0.5) • • • (— 5, 10, — 0.5) in units of A. Magnetic field lines 
are reconnected around the central low-density region. Contrary to the 2D reconnection in 
the x-z plane, reconnection structure involves the x-y plane, because the current sheet is 
highly folded by the RDKI. Plasma inflows mainly come from the indirections, and then 
outflows flow into the indirection. Inflows from the iz-direction are difficult to identify. 
Interestingly, the central reconnection region contains multiple field reversals between the 
folded current sheets. Since plasmas are expelled away from the reconnection region, there 
is a low-density hole around the central region. The reason reconnection takes place in the 
center may be the initial trigger impact. Around the central region, the RDKI is invoked 
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earlier than other locations, and then the folded structure first appears there. Reconnection 
occurs in such a well-developed folded region. In this simulation, the central reconnection 
region is finally overwhelmed by the current sheet corruption, and then the system evolves 
into the turbulent state at t/r c = 140. This is consistent with the 2D RDKI picture. The x-y 
cross sections of the plasma density structure in the upper wall indicate the plasma mixing 
across the double periodic boundaries. Eventually, most energy in the system is converted 
into plasma heat, 2 times hotter than the original state. Particle acceleration by the RDKI 
or reconnection is negligible. 

5. 3D EVOLUTION WITH GUIDE FIELD 

Next, we study a 3D evolution of the current sheet with a guide field in run 3D-B. The 
guide field amplitude is \B y /B \ = 0.5. All other conditions are the same as those of run 
3D-A in £j4j Throughout the simulation run (0 < t/r c < 220), the total energy is conserved 
within an error of 0.6%. 



5.1. Linear Evolution: Relativistic Drift Sausage/Kink Tearing Instability 

First, we study the linear structure of the mode in detail. Compared with run 3D- A, the 
current sheet seems to be more stable in the early stage. It takes t/r c ~ 100 until we observe 
visible changes in the current sheet. At tjr c = 120, we observe a purely growing flute like 
mode in the oblique direction on the upper side of the current sheet in Figure [7]a. We think 
this is a generalized mode between the relativistic tearing instability and the RDKI-type 
instabilities (ZH05b). The mode is (1,1) or its wavevector is fc x A = (0.25,0.25). The two 
panels in Figure [8] show the plasma density slices at z = ±A. We see the structure with 
ki on the upper side of the current sheet at z — A, while we observe another oblique mode 
(1, —1) : k 2 \ = (0.25, —0.25) on the lower side of the current sheet at z = —A. On both 
sides, it seems that the oblique lines are slightly disconnected around (x,y) ~ (±12.8A,0). 
Probably the fastest growing modes do not match the system size, and then they adjust into 
the nearest periodic mode (1,±1). Therefore, it takes a longer time before we observe the 
visible change. The most powerful Fourier modes are (1,±1) in the simulation data, and 
Figure [9] shows ^-profiles of the density perturbation of the two modes. The amplitude of 
the two modes is nearly same, but their peaks are in the other side of the current sheet. We 
confirmed that these oblique perturbations are purely growing unstable modes. Carrying 
out a supplemental 2D PIC simulation under the same conditions in that particular angle, 
we observe similar a similar asymmetric profile, mode (2, ±2), while we obtain mode (1, ±1) 
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in run 3D-B. 

In order to study the nature of the instability, we have also calculated eigen profiles 
of the oblique instability in the current sheet by using the relativistic two-fluid theory. We 
assume the perturbation of 5f oc f(z) exp(ik x x + ik y y + u>it) as introduced in §3.1, and 



then we consider the oblique case of k x ,k y ^ 0. In Figure 10, the eigen functions for 
k x X = 0.25, k y \ = 0.25 and B y /B = —0.5 are presented as functions of z. Perturbed 
magnetic field profiles (5B X , 6B y , i5B z ) and electric field profiles (SE X , SE y , i5E z ) employ the 
same normalization. The other three panels display the current profiles (5J x ,5J y ,iSJ z ), the 
bulk velocity perturbations (i5V +x , iSV +y , SV +Z , where = 5v p + 5v e ), and density profiles 
(5D + = 5d p +5d e and 5D_ = 5d p —5d e ) are normalized by their maximum values. The current 
profiles (SJ x ,5J y ,i5J z ) are calculated from other perturbed values: 5D + ,5V_, where 5V- = 
Sv p — Sv e . These profiles are consistent with the simulation data. For example, the density 



perturbation 5D + in Figure 10 is in excellent agreement with the relevant perturbation in 
simulation data in Figure |9j In general, the perturbation profiles are not symmetric nor 
anti-symmetric with z. As observed in 5B X or 5D + , the mode seems to be substantially 
localized in the upper half of z > 0. The other mode (k 2 ; k x X = 0.25, k y X = —0.25) has an 
opposite structure; it is localized in the lower half of z < 0. 



Figure 11 schematically illustrates the structure of the oblique mode (fci). Although we 
called this mode relativistic drift-kink tearing instability (RDKTI) in ZH05b, we find that 
the structure is rather similar to that of the relativistic drift sausage instability (RDSI), a 
cousin mode of the RDKI. Therefore, we shall rename it relativistic drift-sausage tearing 
instability (RDSTI). For simplicity, we assume that the perturbed structure is localized in 



the upper half of the current sheet in the top panel in Figure 11 In addition, the two 
conventional instabilities in antiparallel configuration; the RDSI/RDKI and the relativistic 
tearing instability are illustrated in the bottom panels of Figure [TTJ The small arrows 
represent the perturbed vectors. The phase of the linear perturbation is set to zero at 
the origin. Therefore, the dashed line in the neutral sheet is phase 7r and the imaginary 
perturbations (phases 7r/2 and 3ti/2) are found between the origin and the dashed line. The 
dashed line (phase 7r) in the neutral plane is equivalent to the X-line of the tearing mode. 
If we look at the eigen profiles in the upper vicinity of the neutral plane (z > 0), plasma 
bulk inflow into the X line (— 5v z < 0), the diverging outflows in i5v x ^ y , the perpendicular 
magnetic field i8B z , and the current enhancement J y are all consistent with the those of the 
tearing mode. In addition, roughly speaking, the density SD + is positive in the O-like region 
and negative in the X-like region. The reconnection electric field E y is positive around the 
X line, although it not strong enough to penetrate into the lower half domain. There is also 
E z structure in i5E z , which is between the X- and O-type regions. The E z structure and the 
relevant charge separation in i5D_ are signatures of the guide field tearing mode or guide 
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field reconnection, as presented in Figures |2]c and |2]e. 

Next, we compare the perturbed structure with that of the RDSI/RDKI in the y-z plane. 
We compare the eigen profiles with those of the RDKI and the RDSI, which are presented 
in Figures 14 and 21 in ZH07. In both cases of the RDSTI and the RDSI, the two-peak 
structure in density profile 5D + stands for sausage-type modulation, while reverse peaks 
stand for kink-type modulation. The charge separation and E z structure of the RDSTI is 
well connected to the those of the RDSI. One specific feature of the RDSTI is the perturbed 
plasma flow structure. The flow direction changes as a function of z, and then the flow is 
parallel to the background magnetic field in the topmost layer. Overall, the structure of 
the RDSTI is complicated, but it is well consistent with the relevant 2D instabilities. The 
RDKTI has almost same structure as the RDSTI, except that perturbation is kink-like in 
the lower side of the current sheet in the y-z plane. 

We investigate the eigen modes over various parameters: the guide field amplitude 
B y /B = 0,-0.5,-1,-1.5 and the wavevectors of (0 < k x \ < 1,0 < k y \ < 1). Their 



fastest growth rates (r c o;j) are presented in contour maps in top panels in Figure 12 We 
take Ak Xt y\ = 0.05 so that 21 2 parameters per map are presented. All of the modes are 
purely growing. Due to the mathematical symmetry, we can obtain eigen modes and growth 
rates for (k x , k y ) from their counterparts in i\k x \, \k y \). We repeat the validity of our linear 
theory again. As discussed in §3.1, our two-fluid approximation will be valid only in the 
long-wavelength range (e.g. \k y \\ < 0.7 for the RDKI/RDSI). Therefore, the results in 



Figure 12 will be reliable only when \k\\ < 0.7. The bottom panels in Figure 12 show the 
type of the fastest eigen mode. By comparing the peak-structure in their density profiles, 
we classify the obtained fastest eigen modes into the following three types: kink type (A 



in Fig. 12), sausage type (B), and neither of them (C). The tearing mode is classified in 
sausage-type modes, because its density perturbation is symmetric with the neutral plane 
(z = 0). Notice that both kink-type mode and sausage-type mode coexist in the same point, 
and that we discuss the type of the fastest mode in the specific parameter range. It seems 
that the RDKI/RDKTI is faster only in the shorter wavelength region along the fc^-axis. In 
the other region, the RDSI/RDSTI replaces the RDKI/RDKTI. 

In the guide field case of B y /B = —0.5, the RDKI/RDSI along the k y -axis is stabi- 
lized by the magnetic tension of the guide field as discussed in §3. The oblique RDSTI 
modes, which are along the background magnetic field lines, survive instead. Some oblique 
modes grow even faster in the guide field case because they are weakly stabilized by the 
magnetic tension in antiparallel case. As the guide field becomes stronger, the most domi- 
nant modes change direction, in accordance with the background magnetic field lines. Since 
the RDKI/RDSI component of the instability is driven by the fc-aligned component of the 
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current, the growth rates of the oblique modes decrease. On the other hand, growth rates 
of the tearing mode component seems to be rather insensitive to \B y /B \. In the strong 
guide field case of \B y /B \ > 1, we can no longer classify some oblique modes to kink-type 
or sausage-type modes (C in Fig. 12). These modes are highly localized on one side of the 
current sheet, and their perturbation is very small on the other side of the current sheet. 



5.2. Nonlinear Evolution 

After the linear stage, the current sheet becomes very thin at several points, where 
two RDSTI/RDKTI modes compress the current sheet from the upper and lower sides. The 
thinning point evolves into a big plasma hole at the center of the simulation box at t/r c = 170 
in Figure [7^. The secondary magnetic reconnection takes place there. Once reconnection 



breaks up, it continues to grow. Figure 13 shows a 2D slice at y = at t/r c = 170. 



The upper panel shows the typical reconnection structure. The outflow velocity is up to 
0.6c, which is between the antiparallel case (0.8c) and the strong guide field case (0.4c) of 
By/ Bq = —1.5. The typical inflow velocity is ~ 0.1c. The current sheet still looks similar to 
a Sweet-Parker current sheet, but careful observation shows several signatures of the guide 
field reconnection. The bottom panel shows the positron current structure in color and the 
E z structure in contour, respectively. The positron current layer is slightly inclined in a 
clockwise direction. However, the electron current layer is inclined in a counter clockwise 
direction. So, the two current layers coexist in a Sweet-Parker-like current sheet. In addition, 
the vertical electric field E z is not negligible. Its amplitude is ~ 0AB in the right side and 
~ — 0.4B in the left side, while the typical reconnection electric field is E y ~ 0.17-Bo- The 
x component of the inflow electric field is E x ~ ±0.05I?o- These signatures are almost same 
in the relevant 2D run: run 3D-A. The width (in y) of the reconnection region seems to be 
limited by the scale of the system or the scale length of the oblique modes. Figure [7]c shows 
a snapshot at t/r c = 200. The oblique bridges are blown away from the neutral sheet due 
to the intense plasma pressure, and then the dense points around (x, y) ~ (0, ±12. 8A) are 
no longer observed. In addition, plasmas are drawn into the central reconnecting point, and 
there is a plasma hole along the X line: x = 0, z = 0. The typical speed of reconnection jets 
is up to ~ 0.71c. Figure [7]ci is the last snapshot of our simulation at t/r c = 220. The system 
structure is highly turbulent, but we observe the filament-like structures in the ^-directions. 
Considering that reconnection dissipates the field energy of antiparallel magnetic fields, it is 
reasonable that we observe structures which are threaded by the guide field B y . 

After magnetic reconnection occurs, the magnetic energy tends to be converted to the 
nonthermal part of the plasma kinetic energy, due to the particle acceleration around the 
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reconnection region. Figure [14] compares the energy spectra in the system for two 3D runs 



and the relevant 2D runs of guide field reconnection. The initial state of two 3D runs is almost 



similar to the spectrum of run 3D-B at t/r c = 140 (Fig. 14, dotted line). In the case of run 
3D-A (bold line), plasma energy is converted into plasma heat, due to magnetic dissipation 
by the RDKI. The nonthermal tail of the spectrum (e > 30mc 2 ) is slightly enhanced due 
to the particle acceleration. The footpoint of the nonthermal tail is approximately same as 
that of the 2D RDKI case, but the tail itself is not as apparent as in the 2D case. Due to 
the irregularity along the x-direction or other 3D effects, particle acceleration by the RDKI 
works less effectively than the ideal 2D case. In the case of run 3D-B, the energy spectrum 
is almost unchanged until relativistic reconnection breaks up. At t/r c = 140 (Fig. [l4j dotted 
line), plasmas are slightly heated by 3% from the initial state. In the late stage of t/r c = 220 



(Fig. 14, dash-dotted line), the nonthermal tail is enhanced due to the particle acceleration 
by reconnection. The nonthermal slope is well-described by the power law with the index 
of —2.8 in a range of 8 < e/mc 2 < 20. This is nearly same as the spectral index in the 2D 
runs. The power-law spectral index is ~ —2.7 at t/r c = 200 in run R3-A, and ~ —2.9 at 
t/r c = 200 in run R3-C. 

The time history of the nonthermal ratio parameter, the ratio of plasma nonthermal 



energy to plasma kinetic energy, is presented in Figure 15 In the case of run 3D-A (Fig. 15 
thick line), there is a small peak after t/r c = 80 due to the dc acceleration by the RDKI, 
but eventually the nonthermal ratio is less than 2%. This is consistent with 2D simulations 
on the RDKI, which reports the nonthermal ratio of far less than 5% (ZH07). However, 



in the case of run 3D-B (Fig. 15, thick dashed line), shortly after the central reconnection 
region appears at t/r c = 170, more than 14% of the kinetic energy consists of the nonthermal 
energy. This ratio is approximately half of the value in the relevant 2D run (run R3-A; Fig. 



15 thin gray line). We think this is due to the limited volume of reconnection region in 
3D configuration. For example, along the neutral line (X = 0) at t/r c = 200 (Fig. uc), 
reconnection and the relevant particle acceleration is active around the center y ~ 0, while 
we do not observe reconnection flow structure around y/X ~ ±12.8. Since reconnection can 
utilize half of the system volume along the y direction, the nonthermal ratio is small in run 
3D-B, compared with the 2D counterpart (run R3-A). If we employ larger simulation box to 
reduce the periodic limitation in the ^-direction, we may observe wider reconnection region 
and then the nonthermal ratio may increase. Anyway, run 3D-B more efficiently generates 
the nonthermal energy than the antiparallel case of run 3D-A, because underlying physical 
mechanism is completely different. Most of magnetic energy is dissipated into the thermal 
energy by the RDKI in the anti-parallel case (run 3D-A). However, in the guide field case (run 
3D-B) a substantial amount of magnetic energy is dissipated into the nonthermal component 
of plasma kinetic energy, associated with magnetic reconnection. 
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6. DISCUSSION 

First, we discuss the system evolution in an antiparallel configuration. We are interested 
in which process dominates, the reconnection or the RDKI, because it greatly changes the 
energy distribution in the system; reconnection involves nonthermal particle acceleration, 
but the RDKI leads to plasma heating. Comparison of linear growth rates (ZH07) suggested 
that the RDKI dominates in the relativistic regime of T/mc 2 > 1. Our results in §4 basically 
support this argument. Although we imposed the external trigger force, it was not strong 
enough to evoke reconnection before the RDKI modulates the current sheet. The dc accelera- 
tion by relativistic reconnection is not likely to evolve; then magnetic dissipation and plasma 
heating by the RDKI would be the main signature of a relativistic current sheet. However, 
we discovered small reconnection regions inside the folded current sheet structure. The re- 
connection generates a density hole around the center, and then the density hall may evolve 
into a global reconnection structure in the larger system. It is true that the RDKI grows 
faster than the reconnection in our relativistic regime of T/mc 2 ~ 1, but once reconnection 
is initiated, we do not know whether it overwhelms the outside RDKI structures. 

Next, we discuss the guide field effect on the system evolution. We showed that the 
RDKI is completely stabilized by the finite amount of the guide field in §3. We do not yet 
understand how much guide field is necessary to stabilize the RDKI, but we expect that 
reconnection mode dominates again under the guide field conditions. In §5, we showed that 
the evolution is more complicated than expected. We considered the oblique instabilities 
(RDSTI/RDKTI) to understand the linear evolution of the current sheet. They grow in 
two oblique directions (e.g. k\ and fc 2 ), which can be interpreted as twin extensions of the 
conventional RDSI/RDKI. In the antiparallel case the kink-type mode dominates, and there 
were no signatures of sausage-type modes. In the guide field case, we found that the sausage- 
type branch (RDSTI) dominates in run 3D-B. In some sense this is quite reasonable, because 
both the RDSI and the tearing instability modulate the current sheets in a symmetric way 
with the neutral plane, while the RDKI modulates the current sheets in an asymmetrical 
way. On the other hand, since the RDSTI/RDKTI have asymmetric structure and since 
they always appear as twins in the guide field cases, it does not matter whether the oblique 
modes are sausage-like or kink-like. The important point is that the twin oblique modes, 
whose wave fronts are parallel to the twisted background magnetic fields, initiates magnetic 
reconnection. Since relativistic guide field reconnection involves particle acceleration, the 
guide field turns on the nonthermal particle acceleration in the 3D system. When we impose 
a stronger guide field, one can see that an angle between two RDKTI/RDSTI branches (k\ 
and k 2 ) becomes wider in accordance with the lobe magnetic field lines. The RDSTI/RDKTI 
slow down, but their growth rates are still faster than the relativistic tearing instability (Fig. 



12). Under the extreme guide field condition, since the dominant RDSTI/RDKTI modes 
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are inclined to the ±x-directions, coupling between two oblique waves may directly lead to 
quasi-2D growth of the tearing instability. 

On the system size limitation, the late-time evolution of our 3D runs may be somewhat 
artificial, due to the periodic boundary effects. In run 3D-A, magnetic diffusion is enhanced 
by the plasma transport across the double periodic boundaries in z. In run 3D-B, the oblique 
mode (RDSTI) seems to be bounded by the system length. Since we see the remnant of the 
other modes (Fig. [8]), we probably failed to observe the most unstable modes. However, as 
long as similar oblique modes dominate they will trigger the secondary magnetic reconnection 
in the similar way. If we set larger simulation box, do we observe multiple reconnection points 
in accordance with the spatial structure of the RDSTI/RDKTI? If so, how are particles 
accelerated in multiple reconnection regions? Or does the X line extend and then does the 
system evolve into the 2D reconnection? Furthermore, recent and ongoing works on large- 
scale evolution of magnetic reconnection exhibit more dynamical behaviors than expected 
(e.g. plasmoid formation and collisions; Daughton & Karimabadi ( 2007[ )). Large-scale, 3D 



evolution of relativistic pair plasma reconnection still remains an open issue, regardless of 
the presence of the guide field. Future simulations may reveal various dynamical behaviours 
beyond our linear and early nonlinear results. 

Regarding the particle acceleration in the guide field reconnection, we studied particle 
acceleration by the parallel electric field E\\ in the 2D case, and we observed similar accel- 
eration signatures in the 3D case, too. In general, as long as we investigated we observed 
nonthermal energy spectra with indexes of —2 ~ —3 by magnetic reconnection, regardless 
of the amplitude of guide field, in both 2D and 3D cases. The upper limit of the acceler- 
ated energy remains to be solved, or it seems unlimited as long as reconnection continues 
to consume magnetic energy. We confirmed that maximum lepton energy easily exceeds the 
Lorentz factor of 80-200 (40-100MeV; see Table [I]). This will also be influenced by large-scale 
evolution of the system. 

Let us briefly discuss potential astrophysical applications, based on physical insights 



from our results. In pulsar winds with relativistically hot pair plasmas (Coroniti 1990 



Lyubarsky & Kirk 2001; Kirk & Skjaeraasen 2003), it is unlikely that current sheets contain 



substantial amount of guide fields because magnetic field lines are highly striped or toroidal. 
Therefore, plasma heating by the RDKI is the most likely process, unless plasma temper- 
ature drops down to nonrelativistic one. On the contrary, in the AGN context, magnetic 
reconnection is quite likely to involve the guide field component. Several authors demon- 
strated electron acceleration due to the field-aligned electric field E\\ by means of test particle 
simulations ( Schopper et al.||1998 Nodes et al.||2003 ). We demonstrated that relativistic par- 



ticle acceleration by E\\ in a fully self-consistent way, including feedbacks from accelerated 
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electrons and other kinetic effects, and so our results provide a theoretical proof of the ultra 
relativistic particle acceleration in the MeV/GeV range, at least in electron-positron pair 
plasma reconnection. In the case of soft gamma repeaters, it is unclear how giant flares 



occur at magnetars. The present models discuss the crustal "quakes" (Thompson & Duncan 



1995, 2001) or the flux tube corruption in analogy with solar flares and coronal mass ejec- 



tions to trigger giant flares. In addition, the tearing instability in the force-free magnetar 



corona is proposed to explain the subsequent bursting activities (Lyutikov 2003). In these 
models, magnetic energy is stored in a magnetic spiral in the star core or the flux lopes, or 
coronal magnetic shear. Therefore, once a flare occurs, it is likely that magnetic reconnection 
involves a magnetic shear or out-of-plane magnetic field. The guide field reconnection will 
occur, and reconnection will be a yet another favorable source of nonthermal particles as 
well as ultra-relativistic shock fronts. 

Recently, relativistic MHD models have been developed to investigate astrophysical 
plasma problems. Our results warn that an MHD approximation is no longer valid in the 
case of relativistic guide field reconnection. As plasma outflow becomes an order of c, charge 
separation structure becomes apparent in the outflow region, and then it breaks down the 
charge neutral assumption of one-fluid theory. Instead, multi-fluid simulations, which deal 
with positively charged fluids and negatively charged fluids independently, are favorable to 
study with guide field reconnection problems. 

Let us summarize this paper. First, we investigated how the guide field affects 2D 
instabilities; the RDKI is stabilized by a finite amount of the guide field, while the recon- 
nection/tearing mode is rather insensitive to the guide field. Then, we studied the nonlinear 
evolution of relativistic guide field reconnection. Characteristic field structure and particle 
acceleration process were investigated. Next, we studied 3D evolution of the current sheet. 
As predicted by 2D studies, the RDKI dominates and dissipates the magnetic energy, but we 
also discovered that reconnection occurs inside the folded current sheet. Finally, we studied 
3D evolution with a guide field condition. The properties of the oblique RDSTI/RDKTI 
mode and the nonlinear evolution of secondary reconnection is discussed, in association with 
the 2D counterparts. Due to the guide field, nonthermal particle acceleration, which gen- 
erates power-law energy spectra with an index of —2 ~ —3, occurs in the relativistic pair 
plasmas. Our results show that the guide field reconnection is a favorable acceleration engine 
in high-energy astrophysical plasmas. 
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Guide Field (IBy/Bol) 

Fig. 1. — Guide field dependence (\B y /B \) of the growth rate u>i, normalized in the light 
transit time r c . The fastest eigen growth rate of the relativistic tearing instability, linear 
growth rate of the fastest growing modes in relativistic magnetic reconnection, eigen growth 
rate of typical RDKI, and linear growth rate of the RDKI are shown. 
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Fig. 2. — Snapshots of run R3-C (\B y /B \ = 1.5) at t/r c = 100.0. (a) Plasma density and 
flows, (b) positron current density J yp , (c) charge non-neutrality [a = [n p — n e ] / [n p + n e }}, (d) 
reconnection electric field E y , and (e) vertical electric field E z are presented, respectively. 
The background counter lines show the magnetic field lines in the x-z plane. 
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Fig. 4.— (a) Spatial distribution of high-energy particles (gray points; e > AOmc 2 ) at 
t/r c = 100. Two typical trajectories are projected in the same plane. The signs show the 
relevant positions at t/r c = 100. (b) Trajectories in the x-y plane. 
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Fig. 5. — Snapshots of the current sheet in run 3D-A at (a) t/r c = 80, (b) 110 and (c) 
140. The gray surfaces show the density surface of (a) n = 2/3n , (b) n = l/3n , and 
(c) n = l/5n , respectively. The plasma density at the neutral plane (z = 0) is projected 
into the bottom wall, with color from black (empty) through blue (sparse) to red (dense; 
n ~ 1.2no)- Panel (d) is a zoomed-in view around the neutral plane (—3.2 < z/X < 3.2) at 
t/r c = 110. The gray lines trace the magnetic field lines. The plasma flow in the 2D planes 
(x = 0,y = 0, and z — 0) are presented as white arrows in the three walls. The light speed 
(v = c) is projected to the arrow length 4. 
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Fig. 6. — The E y structure in run 3D-A at t/r c = 80. 
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Fig. 7.— Snapshots of the current sheet in run 3D-B with a guide field configuration 
(B y /B = -0.5) at (a) t/r c = 120, (b) 170, (c) 200 and (d) 220. Gray surfaces show the 
plasma density of (a) n = 2/3n , (b) 2/3n , (c) l/2n and (d) l/3no, respectively. The 
bottom walls show the plasma density structure in/under the neutral plane at (a-b) z = —X 
or at (c-d) z — 0. 
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Fig. 8. — Slices of the simulation domain in run 3D-B at t/r c = 120. The plasma density is 
represented by gray shading. 
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Fig. 9.— The ^-profiles of the density perturbation in run 3D-B at t/r c = 120. Two 
dominant modes (1,1) and (1,-1) are presented. 
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Fig. 10. — Typical eigen functions of the relativistic drift-sausage tearing instability as a 
function of z, for the normalized wavenumbers k x X = 0.25 and k y X = 0.25. Perturbed mag- 
netic fields: SB x ,SB y ,iSB z ; electric fields: 5E x ,5E y ,i5E z ; electric currents: 5J x ,5J y ,iSJ z ; 
bulk velocities: iSV +x ,iSV +y ,SV +z ; and density fluctuations: 5D± = 5d p ±5d e are presented, 
respectively. 
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3D Relativistic Drift-Sausage (or Drift-Kink) Tearing Instability 
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Fig. 11.— Schematic illustration of the 3D relativistic drift-sausage tearing instability. 
Bottom panels show the relevance to the conventional 2D instabilities. 
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Fig. 12. — Top: Growth rate (r c c<jj) of the instabilities in wavevector spaces of (0 < k x \ < 
1, < k y \ < 1), as a function of the guide field amplitude B y /B = a = 0, —0.5, —1, —1.5. 
Bottom: Three classes of the unstable modes are mapped; A: kink-type modes; B: sausage- 
type modes; and C: others (neither or intermediate). 




Fig. 13. — Snapshots of 2D cut plane at y — in run 3D-B at t/r c = 170. (a) Plasma 
density, flows, 2D magnetic field lines and (6) positron current J py (colored shading), vertical 
electric field E z (contours) are presented. The contour step AE Z is set to 0.1B . 
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Fig. 14. — Energy spectra of 3D/2D runs. Plasma count number for 3D runs are re- arranged 
for comparison with 2D runs. 




Fig. 15. — Time histories of nonthermal ratio parameters for reconnection runs; runs R3, 
R3-A, R3-C, 3D-A and 3D-B. 
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